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SUMMARY 

The purpose of this study was the examination of an empiri- 
cal method of obtaining aircraft trimsfer functions directly from 
in-flight records without recourse to lengthy frequency domain 
techniques. 

The proposed method utilizes a precision analog computer 
and plotter to generate response curves matching the in-flight 
records. 

The computer set-up used for simulation is so arranged 
that each numerator coefficient of the transfer function is directly 
related to a potentiometer setting. At the most, two potentiometer 
settings are unknown, since the first and last coefficients can be 
obtained from the initial slope and steady state value by using 
Laplace initial and final value theory. In the case of longitudinal 
motion the period of the phugoid is of such a long duration as to 
make it impractical to obtain final value steady state conditions 
in actual flight records. This is not prohibitive for this problem 
solution method, however, since one coefficient can be determined 
from the flight record initial slope or suitable accelerometer 
measurement and then the two remaining coefficients are de- 
termined by adjusting potentiometer settings. 

Utilization of this method indicated that errors would be 
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of the order of five per cent in magnitude with the combination 



of computer and recorder used. It is felt that with a high speed 
repetitive computer and a large high-persistence cathode-ray 
oscilloscope, the time of solution and degree of error might be 
reduced if sufficiently accurate flight data were obtained 
initially. 

Further study is indicated for inputs differing from step 
functions such as input with the Laplace transform; 



The theory for tliis case is mentioned briefly in the Proposed 
Technique section. 
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INTRODUCTION 



Considerable theoretical work has been done with frequency 
domain techniques as applied to aircraft dynamic analysis, whereas 
little effort is indicated in the field of transfer function determination 
directly from in-flight transient response data. The purpose of this 
thesis is to investigate a simplified method of determining transfer 
functions by the latter means. 

The empirical transfer function solutions are checked 
against the solutions of the equations of motion utilizing stability 
derivatives and other data of the F-86D aircraft. This aircraft 
was chosen since it is of relatively recent vintage and has been 
in existence long enough to supply reliable information sufficient 
to verify the techniques developed in this thesis. Transfer 
function poles are readily determined from flight test records; 
however, the zeroes require tedious computation when frequency 
domain methods are used. The values of the transfer function 
poles were taken from analytical computations, however, they 
were checked on the computer for validity and accuracy. The 
previous precaution was taken in order to narrow the field of 
investigation and source of error in determining the zeroes which 
were the crux of the overall problem. It was subsequently de- 
termined that this precaution was not necessary and would not 



be required in an actual practical situation. 
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Two basic theorems of the Laplace transformation (See 
Ref. 1) constitute an important part of the development; 

(a) Final value theorem - If the function f(t) and its first 
derivative are Laplace ti-ansformable and if 

L[f(t)] = F(s) 

then 

Limit sF(s) = limit f(t) 
s ->“0 t — r* exo 

(b) Initial value theorem - If the funetion f(t) and its first 
derivative are Laplaee transformable and if 

L[f(t)] = F(s) 

then 

Limit sF(s) limit f(t) 
s — ►©© t 0 

The applieation of these theorems will be shown in eon- 

junction with measurements taken from the assumed in-flight 

transient response curves. 

This study was undertaken by R. W, Ulgen and W. P. 
Vosscler under the guidance aiid direction of Professor R. M. 
Howe as a joint thesis project to satisfy, in part, the require- 
ments for a Master of Science Degree in Aeronautical Engineering 
at the University of Michigan, Ann Arbor, Michigan, during the 



period September 1957 - June 1958. 
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SYMBOLS AND NOTATION 



CX 

b 



c 
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mac 

N 



dCj^/ do^ 

Se 
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L 



m 

M 
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s 



V 

w 



/3 



n 

r 



Angle of attack, degi’ces or radians 
Wing span, feet or inches 
Chord length, feet or inches 
Drag coefficient, D/qS 
Lift coefficient, L/qS 

Moment coefficient about aerodynamic center 

Normal force coefficient, N/qS 

Normal force coefficient curve slope due to 
angle of attack dCj^/dc< 

Pitch moment coefficient, m/qSd 

Slope of the lift curve 

Elevator deflection, degrees or radians 

2 

Acceleration of gravity, ft/sec 

Altitude, feet 

Lift, pounds 

Mass of airplane, slugs 

Mach number 

Dynamic pressure, ^ V /2 
Wing area, square feet 
Airspeed, velocity, knots, mph, ft/sec 
Weight of aircraft, pounds 

Angle of attack in yaw plane, degrees or radians 
Natural frequency, rad/ sec 
Damping coefficient 
Time constant, see 



AIRCRAFT AXIS SYSTEM 







Axis 
Force 
Moment 
Linear Velocity 
Linear Acceleration 
Angular Displacement 
Angular Velocity 
Angular Acceleration 
Moment of Inertia 
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Vp = iU+ jV-»- kvv space linear velocity vector 

Si = iP-t-jQ+kR space angular velocity vector 
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Longitudinal Response Curves 
Longitudinal Simulator Diagram 
Longitudinal Simulator Response Curves 
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PROPOSED TECHNIQUE 



THEORY 

Transfer functions can be obtained from the equations of 
motion which are ratios of poljTiomials of the general form: 



X(s) _ 4. bQ_]^s^ ^ 4- 



m> n 



(f(s) anis^ ajn-ism-1 t ••••-• Uq 
M ultiplying the above by and employing the Laplace 

initial value theorem for a step input: 



Limit s(s”^“”)X = limit 
s >00 t — 



iin- 



nx(t) 



dt 






The slope of the (m-n- 1 ) derivative response curve is then 
equal to b^^ 

The final value theorem for a step input gives: 



limit sX = limit X(t) = — ^ — o 
s — ►0 t — >00 

The first and last coefficients of the numerator can now 
be easily evaluated since the denominator coefficients are found 
by the simple means of measuring damping and frequency. For 
n = 2 there are (n-1) unknown coefficients remaining. 

If the above transfer is mechanized (n- 1 ) potentiometers 
representing the unknown coefficients can be varied so as to 
obtain results matching the in-flight records. This is not overly 
difficult to accomplish since there are a maximum of two 



potentiometer settings which are unknown and these are known 
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at least to an order of magnitude. Tliis would apply to rigid- 
body aircraft. 

A suggested method of procedure for a step function in- 
put is to match phase fairly closely and then match initial slopes 
for the same steady state value. Iteration brings relatively 
rapid convergence. 

For the case where the input transform S (s) is ^ 

s(sT 1 ) 

the final value theorem gives results identical to the previous case. 
The initial value theorem, however, indicates that 



limit 
t— ^ 0 



din-n-l x(t) 



bn ^ 



dtm-n-1 r 

where is the time constant of the input. 

Only examples having step function inputs will be con- 



sidered in this thesis. 
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EQUIPMENT 

The computer used for tlicse studies was a 20 amplifier 
elccti'onic differential analyzer designed and built by the De- 
pai'tment of Aeronautieal Engineering, University of Miehigan. 

A photograph of the eomputer, along witli the Variplotter, is 
shown in Fig. 1. This eomputer had resistor accuracy of . 02%, 
which was more than required for the present study. The 

5 

computing amplifiers had a gain of approximately 10 and were 
manually balanced. Coefficient potentiometers were set to 
an accuracy of 0. 1%. 

The Pace Associates Model HOOD Variplotter used to 
plot the analogue computer outputs is a table-top size, self- 
balancing, potentiometer type of recorder. It records two 
vai’iablcs simultaneously in the form of a rectangular coordinate 
graph. The static accuracy of both the pen and arm is .075 per 
cent. The dynamic accuracy for both is . 2% for plotting speeds 



up to 15 inches per second. 
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EXAMPLES 

Examples of the determination of zeroes for both longi- 
tudinal and lateral functions will be given. 

I. Longitudinal Transfer Function Zeroes 

Since no suitable flight records were available, flight 
data was assumed from solving the aircraft equations on an 
electronic differential analyzer. 

Under conditions of perturbations from steady level 
flight at Mach . 8 and 30, 000 feet the following equations 
were utilized based on the airplane data in Appendix A. 
u = -,00664u + .0497W - 31,9 -32.2 0 

q 

w= -. 0407U - l.OlSw i- 790.3q - 1. 29 6 - 1. 55 6^ 
q = .000475U - . 0119W - .737q -.413<Tg 

e = q 

Scaling the set of equations for use on the analog computer 

= -. 00664U + . 00497(10w) - . 319(100q) - . 322(100 6) 

= -.0407U - . 1018(10w) + 7. 903(100q) - .0129(1006) -.155(10^^ 
= .0476u - .119(10w) -,737(100ci) -. 0413(10«re) 
r lOOq 

These equations w'erc then mechanized as shown in Fig. 2 
(see Appendix B for all remaining figures). The following scaling 
was used: 

1 fps = . 1 volts 



gives 

ii 

w 

100 q 

100 e 



1 rad = . 1 volts 
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1 deg 



. 1 volts 



1 rad/sec = . 1 volts 



The exact transfer function of Yq^ ^ was found from 
this set of equations to be: 




(s^ -h l,762s^ + I0,16s^ . 0195s +. 0311) 



23.6(s^ .983 s -4- . 00855) 



Since there is very little coupling between the short period and 
the phugoid motion it was considered possible that the poles and 
zeroes of the above transfer would be the same as the poles and 
zeroes of the uncoupled transfer functions. 

The approximate formula for the short period transfer 
function as based on Ref. 2 is 



The approximate formula for the phugoid transfer function 
as obtained from the same source is 



The product of the above two expressions does indeed give 
a very close approximation to Y © ^ which makes a lengthy 

O o 



e “ {. 0978 s2 + . 1885 s 1) 



{1.032 s ^ 1) 




(127 s + 1) 



(305.5 s^ + 2.405 s +• 1) 



Newton Raphson approach (see Ref. 3) or other factorization 



technique unnecessary. 



The parameters of the approximate transfer function 



are then as follows: 
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3. 22 radians cO = , 0572 radiaxis 

np 

1.862 seconds Tp = 109.5 seconds 

,304 V = -0685 

These values were used to assure that the response 

curves obtamed from the mechanization of these cuiwes were 

indeed correct. The computer response curves of Fig. 3 had 

the following parameters 

60 = 3. 22 radians - ,0572 radians 

ns np 

Tg -2 seconds Tp = 112 seconds 

ys = -3 - -07 

The above parameters give the poles immediately in the 
canonieal form. 

The detei'mlnation of zeroes for the transfer function 

Y 0 r was considered next. 

® e 

In the interest of a simpler circuit, the constant term 

(. 0311) in the denominator was assumed to be negiigible. An 

's' was then factored out of the remaining terms giving the transfer 

function of Y r 
q <J e 

/ _ -23. G(s^ T .983s T , 00855) _ 

® e " (s^ + 1.762 s 2 + 10.16s + .0195) 

which has the dimensions of degree per degree or radian per 

radian. 

The denominator is now a cubic with coefficients which 
can be equated to potentiometer settings and amplifier gains as 
follows for the longitudinal simulator computer circuit of Fig. 4: 




!Ts = 
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z 1 . = 10, ^ = .1762 

(X 2 K 1 K 2 K 3 = 10.16, K 3 = 2, 0^2 = .508 

0 < 3 K]^K 2 K 3 K 7 = .0915, (ocK)^ = . 1, CX 3 z .010 

Inasmuch as the steady state of the phugoid is difficult 
to obtain in flight tests, it was decided that a more practical 
method of procedure would be to determine the potentiometer 
setting affecting the initial slope by utilizing flight records and 
to vary the remaining two potentiometers. 

The initial slope potentiometer setting corresponding to 
the coefficient of the s term in the numerator was found in the 
following manner. 

The initial slope is found from Fig. 3a or 3b. Recalling 

that 

V<T = ’^2/^3 

we have from Fig. 3 

b 2 = . 4 rad/sec /degree. 

This was used to obtain the potentiometer setting as 
shown in Appendix C. 

As has been indicated before, steady state for the phugoid 
takes too long to obtain and a cut and try approach is more 
expedient. This leaves potentiometers five and six which must 
be adjusted in the computer setup of Fig. 4. Various combinations 
over the full range of each potentiometer were tried. At each 
setting an output trace of q for a 1 ° elevator step was compared 
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witli the assumed flight record. Notuig the effects and trends of 
these combined settings it was possible to come within 5% of 
what the actual required sotUngs should be after only about 18 
ti'ials. Licreascd magnification using si^ecialized high-persistence 
cathode ray oscilloscopes should increase aecuracy of sueh 

measurements. 

Fig. 5 illustrates the technique in a general way uti- 
lizing the following sequenee of settings; 

Trial Potentiometer 

5 6 



A 


. 5 


.5 


B 


.75 


.5 


C 


1. 0 


.5 


D 


1. 0 


.0 



It is to bo noted that the initial slope remains unaffected 
as might be expected. The latter condition suggests that the 
potentiometer 4 setting could also have been obtained by trial 
and error matching instead of from initial sloi^e measurement. 
Variation of potentiometer 5 affects the time and magnitude of 
the first peak as well as setting the continued level of the entire 
eurve. Variation of potentiometer 6 has a negligible effect on 
the initial peak, however, it raises or lowers the tail of the 
cuiwe thus uifluencing the final or steady state value. 



15 



II. Lateral Transfer Function Zeroes 

Again, since no suitable flight records were available, 
the flight response curves were obtained by solving the air- 
craft equations on an electronic differential analyzer. 

In a manner similar to the longitudinal case, the air- 
plane data in Appendix A was used to obtain 

V = -. 179v + 31. 8p - 796r + 32. 2 ^ 

p = -. 435v - 3. 6p + 1. 295r -f . 313r - 1. 015 cT + . 164 <T 

a r 

1- = . 0121V - . 0389p - . 388r -H . 0765p - . 0146 <T - . 0838 cT 

a r 

^ r P + . 04r 

Scaling the set of equations for use on the analog 
computer gives: 

V = -. 179v V . 318(100p) - 7. 96(100r) . Z22(im<j>) 

lOOp =-.0435v- 3.6(100p) + 1.295(100r) + . 313(100r) - 101, 5<f^ 

t 16.4 cT 

r 

lOOr = 1. 21v - . 0389(100p) - . 388(100r) t . 0765(100p) - 1. 46 cT 

a 

-8. 38 (Tr 

lOOjz^ = lOOp + . 04(100r) 

These equations were then mechanized as shown in Fig. 6. 
The scaling was as follows: 

1 fps = , 1 volt 
1 rad = . 1 volt 



1 deg = . 1 volt 






I 



I 
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A static check was performed by observing the voltage 
on inverters 1, 6, and 11 for various initial condition voltages. 
The response curves are shown in Fig. 7. 

The important transfer functions for this set of equations 
were found to be as follows; 



r = 138(s^ - 2.45s - 9.86) 

D(s) 

YM = l->051(s^ -H .582s -F 10.1) 
^ ^ D(s) 



Y r 
vdr 



Y 



vda 



Y 



rS'r 




= -62.7(s^-t- 4. 06s -.0296) 
D(s) 



= -42(s^ - .8723 -.325 ) 

D(s) 

= .073(s^-h 4.493^ -t- 1.45s ~t-.752) 
D(s) 

= .Q945(s^ -I- .3218^ -F 4.48s + 4.51) 
D(s) 



where 

D(s) = s^ -t- 4, 16s^ + 10.6s^+ 38,5s -. 01975. 



Since the transfer functions for uncoupled rolling motion, 
spiral motion, and the dutch roll motion promised to be useful in 
finding the poles and zeroes of the above transfer functions these 
modes of motion were investigated. For uncoupled rolling motion 
we assumed that r = cT^ - 0. This gives 
p = -3.6p - 1.015cf^ 
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or 



i’Ua - 



. 382 



s(. 287s + 1) 

For the dutch roll motion we assumed that p z ^ = cT^ = 0. 
This gives 

V = -. 179v - 796r 
r z .0121V - .388r - .0838 cT 



or 



r- ^ . wooo 

-- 



. 0838(s t .179) 



567s t 9.72) 

For spiral motion we assumed that vzp=rz^z 

rT r 0. This gives 
r 

0 = 179v -t- 31. 8p - 796r + 32. 2 ^ 

0 r -4,35v - 360p t 129. 5r 
0 = 1.21V - 3.89p - 38. 8r 
0 r ^ “P " • 04r 

If we assume in addition that v = p = 0, we get a spiral mode 
time constant of 1875 seconds. 

With these as approximations, synthetic division coupled 
with the Newton-Raphson iteration technique for determining roots 
was used to give 

Y(kr = -.138(s - 3.267) (s + 3.022) 

^ D(s) 

= 1.051(8^ -h .582s 4- 10.1) 

^ D(s) 
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-42(s - 1.15) (s -t- .282) 
D(s) 




. 073(s + 4, 1) (s^ -I- ,32s -t- .2) 
D(S) 




.0945(s t .91) (s^ - ,579s t 5) 
D(s) 



where 



D(s) = (s - .000513) (s t 3.95) (s + .2255s + 9.742) 

The parameters associated with the denominator of the 
various uncoupled transfer functions are then as follows: 

Tj. z ,287 seconds - the time constant of the uncoupled rolling 



T = 1875 seconds - the time constant of the spiral mode 
s 



The parameters associated with the denominator of the 

transfer function are found to be; 

T z , 253 seconds 
r 

T„ z 1950 seconds 
z .324 



motion. 



^ U(jj. z • 324 radians per second - dutch roll 
'fdr = -S’S 




These parameters give the poles immediately in the 
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canonical form. 

The parameters of the simulator circuit shown in Fig. 8 
are determined by setting the coefficients equal to the potentiometer 
settings and amplifier gains as follows; 



«,K,K2 = 4.16 


^^1 


= 1. Kg z 10 


«1 = 


.416 






“ 2, 


o <2 z 


.53 


«jKjK2K3K, = 38.4 


% 


= 2, 


3 = 


.96 



since 

D(s) = s(s^+ 4.16s^ + 10.6s +38.4) 

Two transfer functions, and Y* were investigated, 

as being representative. The value of the potentiometer setting 
corresponding to the highest degree 's' term in the numerator was 
found as in the longitudinal case. See Appendix C for details. The 
remaining potentiometers were adjusted giving the plots and variations 
shown in Fig. 9. 



I 
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RESULTS AND CONCLUSIONS 

By determining potentiometer settings corresponding to 
initial slope and also the steady state value, where available, 
it was possible to obtain accuracy in matching the curves such 
that the zeroes were determined with errors of less than 5%. 

In the case of q/J’^ , where two potentiometer settings 
were unknown, only 18 trials were necessary to obtain this 
accuracy. 

See Fig. 5 and Fig. 9 to obseiwe the variation in re- 
sponse curves for variation in numerator potentiometer settings. 
It is concluded that the best general procedure is as follows: 

1) Obtain estimates of poles and zeroes from pre-flight 
data. 

2) Obtain poles from flight records involving step inputs 

and coefficients of first and last terms in numerator 

corresponding to initial slope and steady state, if 

possible. Then compare with previous estimates. 

Limit d"^"”x/dt”^"“ = ^n ^ 
t — 0 

Limit X(t) = bQ<r /a^ 
t — ►exj 

3) Adjust potentiometers to correspond to the known 
parameters. 

4) Set estimates in those potentiometers which must 



be adjusted. 
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5) Make preliminary adjustments to observe which 
potentiometer has more effect in the case where 
two potentiometer settings are unknown, 

G) Adjust the more critical to give ioose matching, 
observing trend of frequency, phase shift, peak, 
and steady state values, 

7) Adjust the remaining potentiometer for a closer 
match while obseiwing trends, 

8) Return to the first potentiometer and repeat ad- 
justment. 

9) Continue iterating for as close a match as desired. 
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APPENDIX A 



AIRCRAFT DATA AND STABIUTY DERIVATIVES 

The following data is given in references 4 and 5 for 
the F-86D aircraft at 30, 000 feet, Mach 0.8, weight 15, 000 lbs. 



S = 288 ft^ 


c = 8. 09 ft b 


= 8625 slug ft^ 


lyy = 29400 slug ft^ 


1^2 = 35, 300 slug ft^ 


I^x- 2,700 slug ft^ 


V = 796 ft/sec 
P 


Cl = -.82/radian 
Q 


<?< = 2. 3° - 0. 04 rad 




^ , 00089 slug/ft^ 


^Yr = 


= 0. 185 




--6.5 

Mq 

= 4.62/radian 


Ct = -0. lO/radian 

^/3 




Ckt = -117 

^/3 




Ct =-.433 

Dp 




C, = . 15 


Cd« - • 


Cnp = --oi 


Cj) - . 018 


Cm =--2 
Nr 


Cm ^ = -. 0185/degree 


C L = - . 0029/ degree 

^ a 


Cl ^ = 0 - 0089/degree 


Cr^ =. 00043/degree 



^ =- 000055/ degree 
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Fig. 2 

Cumputer Circuit for Solving the Longitudinal Flight Equations 
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Lateral Simulation for sr/Si 
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APPENDIX C 



Determination of Potentiometer Settings 

For the longitudmal case we found that 
b 2 = . 4 rad/secVdegree. 
from Fig. 3 a or 3b. 

The scale that we wish the plotter to have is 
1 inch - . 1 rad/ sec^ 

with the pen scale factor set to 10 volts per inch. 

For a one degree step input the initial q should be equal 
to 40 volts, since 

2 

40 volts = . 4 rad/ sec . 

The equation relating the 40 volts to the parameters of 
cui 



the circuit can be determined by reference to the figure below. 

^ 



. [? / I>(s)] (^4 < ( yKj Kvj 5^^ 
[ 2 /]>Cs)] (®<4K^(re<'syK3K7)s^=X^S) 



Limit sX(s) = Qq = (©< 4 ^^ volts 

s >■ 00 

For ^ =2 volts 

© 

= 2 

In brief then for the initial slope 

(K/ag) = (KvAri) (b2/a3> 

and for steady state 

(K/a^) {cXqK'q Sq) = (Kvi/Kri) (bo/a^) 



C-2 



where 

K = denominator multiplying factor 

K^. - the pen scale factor on the plotter (in the above 
example it was ten volts per inch. ) 

Kj.j z the number of radians /second or the number of 
radians/second /second (in the above example 
it was . 1) 
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